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Abstract 

This paper investigates the pricing and hedging of variance swaps under 
a 3/2 volatility model. Explicit pricing and hedging formulas of variance 
swaps are obtained under the benchmark approach, which only requires 
the existence of the numeraire portfolio. The growth optimal portfolio is 
the numeraire portfolio and used as numeraire together with the real world 



>- ■ probability measure as pricing measure. This pricing concept provides 

minimal prices for variance swaps even when an equivalent risk neutral 
(^ . probability measure does not exist. 
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1 Introduction 



A significant source of risk in a financial market is the uncertainty of the 
volatility of equity indices. During financial turmoil, volatility risk is of 
extreme importance to investors and derivative traders. Additionally, due 
to large and frequent shifts in the volatility of various assets in volatile 
periods, there is a growing practical need for appropriate models that al- 
low for the realistic modeling of volatility, the pricing of related financial 
instruments and the hedging volatility risk. In 1993, the Chicago Board 
Options Exchange (CBOE) introduced a volatility index, the VIX, see Fig- 



ure 1.1, based on the implied volatilities of options on the S&P 500 index, 
see Figure 1.2. The volatility of this diversified index has attracted the 
most attention in the current literature on volatility derivatives. As such, 
this paper singles out as the security of interest a similarly diversified equity 
index and studies its volatility and related derivatives. 




Figure 1.1: VIX from 1990 to 2010. 




Figure 1.2: S&P 500 from 1990 to 2010. 



Variance swaps have been traded in over-the-counter markets since the 
collapse of Long Term Capital Management in late 1998. In particular, 
variance swaps on stock indices have been traded actively as a hedge for 
volatility risk. Investors and fund managers alike, have developed an in- 
terest in volatility derivatives since these instruments may substantially 
increase the value of their holdings, even if the equity market index ex- 
periences a major drawdown. At least theoretically, these derivatives can 
provide some protection against severe market downturns. How effective 
such portfolio insurance is, from a macro-economic view point, remains an 
open question. In particular, when a large and increasing number of pen- 
sion funds, insurance companies and other investors rely on this type of 
insurance, it is not clear whether the sellers of variance swaps will be able 
to serve their obligations in a market crash. 



There exists a substantial literature on volatility modeling, and one may 
refer to Cont and Tankov (2004) as one of many references. Several pa- 
pers, which explore stochastic volatility models, have pointed at the seem- 
ingly undesirable property of some models, where the moments of squared 
volatility of higher order than one may become infinite in finite time. Exam- 
ples are given in Andersen and Piterbarg (2007), Lion and Musiela (2007) 
and Glasserman and Kim (2010). Furthermore, there exist various papers 
discussing the general problem of pricing and hedging variance swaps, in- 
cluding Brenner et al. (1993), Griinbuchler and Longstaff (1996), Carr and 
Madan (1999), Chriss and Morokoff (1999), Demeterfi et al. (1999), Brock- 
haus and Long (2000), Matytsin (2000), Javaheri et al. (2002), Swishchuk 
(2004), Howison et al. (2004), Carr et al. (2005), Windcliff et al. (2006), 
Zhang and Zhu (2006), Zhu and Zhang (2007), Sepp (2007), Elliott, Siu 
and Chan (2007) and Carr and Lee (2009). 

Recently, Carr and Sun (2007), as well as, Itkin and Carr (2009) discussed 
in two interesting papers the pricing of variance swaps under the, so called, 
3/2 volatility model. Earlier versions of 3/2 volatility models were studied, 
for instance, in Cox et al. (1980, 1985), Platen (1997), Ahn and Gao (1999), 
Heston (1999), Lewis (2000), Andreasen (2001), Platen (2001) and Spencer 
(2003). 

The current paper studies the case of a model which it calls the real world 
3/2 volatility model. Explicit pricing formulas of variance swaps are ob- 
tained. This model generalizes the minimal market model (MMM), which 
was introduced by Platen (2001) and does not have an equivalent risk neu- 
tral probability measure. The real world pricing concept of the benchmark 
approach aims to identify the minimal prices for variance swaps. 

The paper is organized as follows: Section 2 introduces variance swaps, and 
relates these to the popular log-contract hedging. Section 3 introduces a 
diversified equity index model. Section 4 demonstrates the estimation of 
parameters. Section 5 studies the risk neutral 3/2 volatility model. Section 
6 considers the real world 3/2 volatility model. Section 7 introduces the real 
world pricing concept. The analytical pricing formulas for variance swaps 
are obtained in Section 8. Section 9 discusses hedging issues. Section 10 
shows some numerical results. Finally, Section 11 concludes. 



2 Variance Swaps 



For simplicity, assume throughout the paper that one has zero interest 
rate. A variance swap is a forward contract on annualized variance, which 
is the square of the realized annualized volatility. Let a^q^ denote the 
realized annualized variance of the log-returns of a diversified equity index 



or related futures over the life of the contract, such that 

1 r^ 

<T--=7^j ^Idu, (2.1) 

where (T„ denotes the volatility of the index or futures at time t. Assume 
that one can trade the futures or index at discrete times tj = iA for i G 
{0, 1, ...} with time step size A > 0. The period A between two successive 
potential trading times is typically the length of one day. If Ff. denotes the 
discounted index price at time t,, then the return G^. for the period before 
this time is defined as 

Ft -Ft 

Gu = % '-' (2.2) 

for i G {1,2,...}. In practice, variance swaps are often based on the realized 
variance, evaluated from daily closing prices, of the form 
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where T = tn = J27=ii^i ~ ti-i) = jji and N is the number of trading 
days per year. Hence, variance swaps with idealized payoffs, depending on 
the realized variance, as defined in (2.1), are only approximations to those 
of the actual contracts. As indicated in (2.3), both returns or log-returns 
may be used to define payoffs on realized variance. Furthermore, by (2.1) 
and (2.3) it becomes clear that the basic elements of this type of payoffs 
are some form of squared volatility to be paid at a given time. This is also 
why the paper will focus on the pricing of a payoff that delivers squared 
volatility. 

As is common in most of the literature, (fi, At^ A, Q) denotes some assumed 
risk neutral, filtered probability space. Here Q is denoting the assumed 
risk neutral probability measure. The filtration A = {At)te[o,T] niodels 
the evolution of the market information over time, where At describes the 
information available at time t G [0,T]. Let K^ denote the delivery price 
for realized variance and L the notional amount of the swap in dollars 
per annualized variance point. Then, the payoff of the variance swap at 
expiration time t„ is given by L{al^^ — K^). Intuitively, the buyer of the 
variance swap will receive L dollars for each point by which the realized 
annual variance 5"q^^ has exceeded K^. Under the risk neutral approach, 
the value of the variance swap can be evaluated as the expectation of its 
discounted payoff with respect to the assumed risk neutral measure Q. This 
value is equivalent to the value of a forward contract on future realized 
variance with strike price K^. 

The risk neutral value V of the variance swap at time t = over the period 



[0, tn] is given by the expression 

V = E^ (lial^ - K,)^ 

= LE^ialJ-LK, (2.4) 

where E'^ denotes expectation under the assumed risk neutral measure Q. 
Hence, the valuation of the variance swap relies on calculating the risk neu- 
tral expectation E^la^^^^) of the realized variance. Note by (2.3) and (2.1) 
that this involves the computation of a sum of risk neutral expectations 
of squared volatility. Therefore, the problem can be reduced to evaluating 
the risk neutral expectation of squared volatility, as mentioned earlier. 

To highlight some of the critical issues with the current methodology, recall 
the popular log-contract hedging approach to realized annualized variance, 
as employed for instance in Neuberger (1990), Dupire (1993), Demeterfi et 
al. (1999) and Carr and Lee (2009). In the above literature the payoff a^^ , 
appearing in (2.3), is usually approximated by the expansion 
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where the return G(. << 1 is assumed to be "small". To be precise in 
our notation, 0{Gf.) represents a term that is smaller in absolute value 
than some constant times |GtJ'' when G^. asymptotically vanishes. The 
first two terms on the right-hand side of (2.5) form the profit and loss from 
a dynamic position in the index and a static position in options on the 

index, respectively. For the dynamic component, one holds M^ ' ^ ^ 



futures contracts on the index from day tj_i to day tj. For the static 
component, one holds -^^dK European put options with strikes K below 
the initial index value Ft^. One also holds ^^dK European call options 
having strikes greater than the initial index value F^g. As such, one holds 
more puts than calls. In the literature it is argued that the third and fourth 
terms on the right hand side of (2.5) represent the most significant source 
of error when approximating the payoff ctq^^. In particular, the fourth 
power of the return appears in the remainder term in (2.5). This paper 
emphasizes that the expansion (2.5) is of a pathwise nature and relies on 
the assumption that |Gi. | is "sufficiently small" for all scenarios, which is 
a rather strong and potentially unrealistic assumption, as we will discuss 
below. 

It was shown in Platen and Rendek (2008), that the Student-t distribu- 
tion with approximately four degrees of freedom is the typical estimated 
log-return distribution of diversified world stock indices when such indices 



are denominated in a currency. Similar empirical evidence is provided in 
Markowitz and Usmen (1996a, 1996b) for the S&P 500. However, it is 
well-known that the fourth moment of the Student-t distribution with four 
degrees of freedom is infinite. In view of the assumed representation (2.5) 
this raises the critical question whether the error truncation in (2.5) is a 
sensible one for practical purposes. It will turn out that the real world 
3/2 volatility model generates Student-t distributed log-returns for the un- 
derlying equity index. In the case of the stylized minimal market model 
(MMM), these log-returns have four degrees of freedom. But also for gen- 
eralizations of the MMM with higher degrees of freedom it is a matter of 
fact that extreme log-returns of diversified equity indices occur from time 
to time, as they arise also in reality rarely but also with some nonvanishing 
likelihood. Under the 3/2 volatility model this happens with a probability 
that is not negligible. Therefore, one may question whether it is prudent 
to assume in (2.5) that \Gt^\ is always "sufficiently small". 

Therefore, when pricing volatility derivatives, one has to deal with expec- 
tations. Hence, the expectation of the final term of (2.5) may create a 
problem in this respect since it can be extremely large or even infinite 
under models that generate Student-t log-returns with realistic degrees of 
freedom. Therefore, based on these facts extreme caution should be taken 
when pricing and hedging volatility derivatives using the type of expansion 
given in (2.5). The authors of the current paper would like to argue that 
the reliance on the assumption that the expectation of the error term in 
(2.5) is "small" should be avoided. Alternative pricing and hedging meth- 
ods should be employed that do not use any such assumption. To visualize 
this problem. Figure (2.3) exhibits the cumulative fourth term on the right 
hand side of (2.5) for S&P 500 index daily log-return data, from 1990 un- 
til 2010. One notes that during the period from September to October 
2008 the cumulative fourth term increased by approximately 300% in just 
a few days. This means that Gf. itself must have experienced extra ordi- 
narily high values. This figure indicates, at least visually, that the fourth 
term in (2.5) cannot be easily neglected as a remainder term that is always 
"sufficiently small". In reality, a severe approximation error has to be ex- 
pected when employing the widely propagated methodology, which uses a 
truncated version of the expansion (2.5). 

To avoid these problems in the valuation of variance swaps, we suggest 
to adopt the benchmark approach, see Platen and Heath (2006). In fact, 
as described in Section 7, the variance swap will be evaluated under this 
approach by using the expectation of its benchmark payoff with respect to 
the real world probability measure P. In principle, this is the value of a 
forward contract on the future realized variance ctq j, with strike price K^. 
Since one considers all quantities under the real world probability measure, 
this excludes the modeling error that potentially arises from the assumption 
on the existence of an equivalent risk neutral probability measure. 




Figure 2.3: The cumulative fourth term, X^ILi ^i^tjy using daily S&P 500 
data from 1990 to 2010. 

3 A Model for a Diversified Equity Index 



Denote by rt the predictable interest rate at time t. The savings account 
value at time t is denoted by S^ and given by the SDE 



dsi 



ndt, 



where r^ may be a stochastic or deterministic interest rate. 

Denote by Ft the value of a discounted diversified equity index at time 
t > 0, and assume that 

Ft = AtYt (3.1) 

for t > 0. Here At stands for the average growth of the index, where 

At = aexp{7]t} (3.2) 

for t > 0, with a > and long term net growth rate t] E ^. The mean 
reverting dynamics of the index around its average exponential increase is 
modeled by an ergodic process Y = {Yt, t > 0}, which is assumed to satisfy 
the SDE 



dYt = (ai - bYt)dt + ^/2^tdW{t) 



(3.3) 



for t > with Yq > 0. Here 6 > is the speed of adjustment, which results 
for Yt in the reference level ^ as long run mean. The process y is a square 



root process of dimension 6 = -^, where we assume 6 > 2. Furthermore, 
W = {W(t),t > 0} is a standard Wiener process. 

By an application of the Ito formula one obtains from (3.1), (3.2) and (3.3) 
for Ft the SDE 



dFt = Ftiixtdt + atdWt) 



(3.4) 



for t > with Fq = AqYq. Here one has the expected rate of return 



fit = r]-b + ^ (3.5) 



and the volatihty 




(^t = \h7- (3-6) 



Therefore, the squared volatihty vt = a^ satisfies under this model the SDE 
dvt = d{ — 



hvt + {l-^)vl]dt-vf^dWt. 



(3.7) 



for t > 0. The dynamics for the volatility of the index turn out to be those 
of a 3/2 volatility model, see Cox et al. (1980, 1985), Platen (1997), Ahn 
and Gao (1999), Heston (1999), Lewis (2000), Andreasen (2001), Platen 
(2001), Spencer (2003), Carr and Sun (2007) and Itkin and Carr (2009). 
Furthermore, they are ergodic, and, thus, due to their stationary density 
remain always within a certain range, even over long periods of time. 

From (3.4) it follows that the market price of risk dt is under the above 
model of the form 



(3.8) 

This shows that the market price of risk becomes large when the market 
index attains relatively small values. As one can see from Figure 1.1 and 
Figure 1.2, this reflects well reality. 



4 Parameters Estimation 



This section shows the estimation of parameters of the model. A set of 
parameters which are needed to be estimated are a, r), ai, b and 7. Figure 
4.4 plots the logarithm of the S&P 500 total return index, observed in US 
dollars from 1920 until March 2010. Furthermore, it also plots a linearly 
regressed function of time: 0.1028t — 1.3929. By a least square flt of the 
logarithm of the S&P500, we flnd a = 0.248 and r] = 0.1028. 

Remarkably, this model reduces to the stylized MMM proposed in Platen 
(2001) when setting parameters ai = 1, b = 1] and 7 = |, see also Platen 
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Figure 4.4: Logarithms of the S&P 500 from 1920 to 2010 as well as linear 
fit. 

and Heath (2006). In this special case the model has only two parameters, 
a and rj. One notes its long term growth with seemingly stationary fluctu- 
ations around the average growth. By taking the logarithm on both sides 
of (3.1) it follows that 



In(Fi) =7]t + \n{Yt) + \n{a). 



(4.1) 



The slope of the fitted line in Figure 4.4 equals the net growth rate r^, which 
turns out to be approximately 10% for the US market when considering the 
last hundred years, see also Dimson, Marsh and Staunton (2002). The value 
In(Y'i) is the logarithm of a stationary process at time t, which is modeled 
in (3.3) as a square root process having a gamma density as stationary 
probability density. This fits well what one observes in Figure 4.4. The 
scaling parameter a can be fitted by ensuring that the average of Yt equals 
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Figure 4.5 shows the normalized S&P500 which is obtained by dividing the 
index by the exponential function given in (3.2) using the above estimated 
values of a and r] respectively. Note the stationary density of the square 
root process in (3.3) has the expression 



p{y) 



r(f) 



exp{--i/} 



(4.2) 



for y G (0, oo), where S > 2, see Platen and Health (2006). 



Figure 4.6 shows the histogram of the normalized S&P 500 together with 
a fit to the above gamma density with estimated dimension 6 = 15. Set 
b = 1] and ai = b, then 7 = 0.0137. 




Figure 4.5: Normalized S&P 500 from 1920 to 2010. 




Figure 4.6: Histogram of the normalized S&P500 with gamma fit. 

5 Risk Neutral 3/2 Volatility Model 



Since the literature on volatility derivatives focus on the risk neutral ap- 
proach we discuss some important results along these lines of research. The 
authors of the two interesting papers: Carr and Sun (2007) and Itkin and 
Carr (2009), derived under some assumed risk neutral probability measure 
a 3/2 volatility model. It is the generic stochastic volatility model that 
naturally emerges under the following two plausible assumptions: the sta- 
tionary volatility ratio hypothesis and the maturity independent diffusion 
hypothesis. The first means that the ratio of the volatility of the variance 
swap rate to the instantaneous volatility of the underlying asset only de- 
pends on the variance swap rate and maturity, but does not depend on 
calendar time. The second assumption requires that the risk-neutral pro- 
cess for the instantaneous variance is a diffusion process whose coefficients 
are independent of the variance swap maturity date. For the detailed de- 
scription of these assumptions and the derivation of the risk neutral 3/2 
volatility model, we refer to Carr and Sun (2007). This model, which will 
be described in more detail further below, appears to reflect well empir- 



ical properties of short term derivatives on realized variance or volatility. 
Some of such empirical evidence is described in Poteshman (1998), Chacko 
and Viceira (1999), Ishida and Engle (2002), Jones (2003), Javaheri (2004), 
Bakshi et al. (2004) and Platen and Rendek (2008). 

Because of the strong empirical support that 3/2 type volatility models 
attract, this paper places this type of model in a broader setting, which 
allows to raise interesting questions from a wider perspective. In particular, 
the dynamics of some 3/2 models is studied under the real world probability 
measure P. For instance, it links the 3/2 volatility model to the stylized 
minimal market model (MMM), described in Platen (2001) and Platen and 
Heath (2006). The stylized MMM represents a 3/2 type volatility model 
when formulated under some assumed risk neutral probability measure, as 
will be shown below. 

As before, denote by Ft the discounted price at time t G [0,T],T < oo, of 
the underlying diversified equity index. The instantaneous variance of this 
index shall be denoted by Vt-, which represents its squared volatility. 

To remain close to the work by Carr and his co-authors, this section mod- 
els the stochastic processes under some assumed risk neutral probability 
measure Q. Assume, as in Carr and Sun (2007), that the equity index Ft, 
or its futures, satisfies the following stochastic differential equation (SDE) 

dFt = Fty^tdZt, (5.1) 

and its squared volatility Vt the SDE 

dvt = lp{t)vt + qvt ]dt + evf^dWt, (5.2) 



t G [0,T]. Here pit), t G [0,T], and q are assumed to be given real valued, 
deterministic quantities. Furthermore, Z = {Zt,t G [0,T]} and W = 
{Wt,t E [0, T]} represent possibly dependent Q-standard Wiener processes. 
The dependence between Z and W is modeled by the covariation 

[Z,W]t = pt, (5.3) 

for t G [0,T], with given constant correlation parameter p G [—1, 1]. Note 
the power 3/2 for the squared volatility Vt in the diffusion coefficient of the 
SDE (3.7) from which the model derives its name. 

It is well-known that the movements of the squared volatility of a diversified 
equity index are in reality negatively correlated to those of the index itself. 
This stylized empirical fact has been called leverage effect, see Black (1976). 
Realistic calibrations identify a significantly negative correlation parameter 
p. Therefore, before going into the analysis of a two-factor model, as studied 
in Carr and Sun (2007), this paper focuses on the key properties of a one- 
factor version. To achieve this, the paper sets p = —1. Thus, a single 



Wiener process is driving the index, as well as, its volatility. The risk 
neutral model leaves some ambiguity about the real world dynamics of 
the diversified index and its volatility. To fix this one can consider the 
one-factor model that follows from (5.1) by setting 

., = ?2^^^^pM. (5.4) 

Ft 

Here a > is a scaling parameter, and rj > plays the role of the long 
term average growth rate of the index Ft under the real world probability 
measure P. As mentioned earlier, the model refiects in a simple manner 
the empirical fact that the volatility increases when the index decreases. 

The relation (5.4) yields, by (5.1), the SDE 

dFt = Ft^^^^0^dZt 
yi't 



V^ 



= y^exp{jt}^tdZt (5.5) 

for t > 0. Under the assumed risk neutral measure Q this is the SDE of a 
time transformed squared Bessel process of dimension zero, see Revuz and 
Yor (1999). One confirms by application of the Ito formula to (5.4) that 
by using (5.5) one obtains for Vt a 3/2 volatility model with SDE 

dvt = di 270; 

V Ft 

= _ 275expM ^^^^^^^^ ^ 27aexp{r^t} ^^^^^^ 27aexp{^t} ^^ 

= Vtir] + vtjdt-v^ dZt (5.6) 

for t G [0,T] with vt > 0. Consequently, the squared volatility satisfies the 
SDE (5.2) with the parameter choice p = —1, p(t) = r], q = 1 and e = 1. 
Of course, there exist other ways of assuming the real world dynamics of 
Vt other than assuming the functional relationship (5.4). However, such 
relationship appears to be reasonably plausible and rather tractable. 

The squared volatility process, characterized by the SDE (5.6), has a non- 
linear drift and a nonlinear diffusion coefficient. At a first glance, this seems 
to make it difficult to obtain statements about its probabilistic properties. 
However, due to the fact that Ft in (5.5) is a squared Bessel process, which 
is an extremely well studied stochastic process, one knows the analytic form 
of the probability density of Vt and many other properties. As in Carr and 
Sun (2007) one may now study the expression 

R. = ^, (5.7) 

Vt 



which equals in our case the value Rt of a square root process with 

Rt = Ftexp{-r]t}. 
More precisely, one obtains by the Ito formula the SDE 



dRt 



2'-fa 



vt[r] + vt]dt- v^ dZt 



+ —^vfdt 



^^^ - 27a + 2-fa)dt + 2-fav;^^^dZt 



Vt 

= -r]Rtdt+ ^2-fdiRtdZt. (5.8) 

It is no surprise that under the assumed risk neutral probability measure 
Q, not only the squared Bessel process F but also the square root process 
R = {Rt,t G [0,00)} is of dimension 6 = 0, see Revuz and Yor (1999). 
Furthermore, according to the just mentioned reference, this process hits 
the level zero with strictly positive Q-probability in any nonzero finite 
time period. In such an event the volatility explodes, see (5.7). This 
gives an indication that under Q the expectation of the squared volatility 
Vt = -^ rnay potentially become infinite. Indeed, under the assumed 
risk neutral probability measure, the first negative moment of the squared 
Bessel process R of dimension 5 = is infinite, see Revuz and Yor (1999). 
Consequently, the risk neutral expectation of squared volatility has under 
the above risk neutral 3/2 volatility model an infinite value, that is, 

E'^ivr) = 00. (5.9) 

Thus, the typical building block of the payoff of a variance swap on a 
diversified equity index does not appear to have a finite risk neutral price 
under the above risk neutral 3/2 volatility model. 

At this point it is important to mention that the assumptions of (5.6) in 
relation to the parameter choice for the risk neutral 3/2 volatility model 
violate those imposed in Carr and Sun (2007) and Itkin and Carr (2009). In 
particular, these authors require in the SDE (5.2) that the relation q < \ 
holds. As a consequence, the dimension 5 of their square root process R 
under the assumed risk neutral measure Q is greater than two. This condi- 
tion guarantees that the volatility will not explode in finite time under the 
assumed risk neutral probability measure Q. The current paper considers 
a different model that is not covered by Carr with his coauthors. It is ar- 
gued that there is no need to restrict the volatility dynamics such that it 
avoids a volatility explosion under some assumed risk neutral probability 
measure. What really matters are the dynamics of the volatility process 
under the real world probability measure P. Under this probability mea- 
sure the volatility should be realistically modeled and, thus, not explode, 
as one observes in reality. As mentioned in the introduction, also other 
volatility models, in particular those that model the leverage effect, may 



create similar volatility explosions under an assumed risk neutral probabil- 
ity measure resulting in infinite prices, see Andersen and Piterbarg (2007), 
Lion and Musiela (2007) and Glasserman and Kim (2010). 

It has been shown in Platen and Heath (2006), and it is also discussed in 
Section 7, that the formally obtained risk neutral price of a nonnegative 
contingent claim can be significantly higher than the real world price to be 
identified in Section 7. This means, if there would be a volatility explosion 
under an assumed risk neutral measure there may still be no volatility 
explosion under the real world probability measure. Thus, there is a chance 
that the real world price of squared volatihty could be finite for a 3/2 
volatility model. 



6 Real World 3/2 Volatility Model 



To study the dynamics of the underlying diversified equity index and its 
volatility under the real world probability measure, this paper adopts the 
benchmark approach described in Platen and Heath (2006). This approach 
does not require the assumption on the existence of an equivalent risk neu- 
tral probability measure. It is only assuming the existence of the numeraire 
portfolio for the given market, see Long (1990), Becherer (2001), Karatzas 
and Kardaras (2007) and Kardaras and Platen (2008). If there exists an 
equivalent risk neutral probability measure for a given model, then the 
benchmark approach recovers fully the results of the risk neutral approach. 
Otherwise, it still provides a sound derivative pricing methodology, as will 
be explained in Section 7. 

Under the real world probability measure P one needs also to model the 
drift in the SDE for the real world dynamics of the diversified index Ft. 
This is not necessary under the risk neutral approach and long term trends 
do not play any role under the classical methodology. 

By (3.4), (3.5) and (3.6) the SDE for Ft under the real world probability 
measure P has the form: 

dFt = Ft\(r]-h + ^)dt + ^tdzX (6.1) 

Here 

Zt = Zt- f 9Js, (6.2) 

Jo 

forms a Wiener process under the real world probability measure P, and 6t, 
see (3.8), is the market price of risk, which provides the link to the process 
Z. The squared volatility for the above 3/2 volatility model satisfies by 



(3.7) and (6.2) the SDE 

dvt = vt(b + {l- ^)v^ dt - vl'''dZt (6.3) 

for t G [0,T], where Z is the Wiener process as given in (6.2) under P. 
The real world dynamics of the process Y = {Yt = — , t > 0} is given by 

r,^?2.^.-« , (6.4) 

Vt a aexpjr^tl 

and satisfies the SDE (3.3). Obviously, the square root process Y has 
a stationary density. Furthermore, it is known that it never hits zero. 
Therefore, the above 3/2 volatility model has no volatility explosion caused 
by Yt hitting zero. Additionally, the parameter rj can be interpreted as the 
long term growth rate of the discounted diversified index, which is a key 
macro-economic variable. 

To elaborate on the interpretation of the real world dynamics of the diver- 
sified index one notes from (3.1), (3.2) and (3.3) that the index Ft satisfies 
the SDE 

dFt = I (^ — b)Ft + aia exp{r]t} ]dt + ^/2'ya exp{rit}FtdZt, (6.5) 

where Z is a P- Wiener process. Essentially, the deterministic drift in (6.7) 
models the increase per unit of time in the underlying "fundamental" value 
of the equity index Ft. On the other hand, the remaining martingale term 
in (6.7) reflects the speculative fluctuations of the index. This parsimonious 
model makes good economic sense, in particular, in the very long term. 

Under the real world probability measure the above model appears to be 
a reasonable model for the dynamics of a well diversified equity index, e.g. 
the discounted S&P 500 total return index. Carr and Sun (2007) and Itkin 
and Carr (2009) use arguments from variance swap modeling to derive in 
a plausible manner a risk neutral 3/2 volatility model. 

Recall that under the real world probability measure P the above square 
root process Y = {Yt,t G [0,T]} has dimension 5 = — > 2, a process 
with stationary density, never reaching zero under P. Under an assumed 
risk neutral probability measure Q the dimension 6 of R, and thus the 
dimension of the square root process Y, is zero. Consequently, there is a 
strictly positive risk neutral probability for the event that this process hits 
zero in a finite time period, see Revuz and Yor (1999). As a result, the 
measures Q and P do not have the same events of measure zero. This 
leads to the conclusion that under the above version of the 3/2 volatility 
model the putative risk neutral measure Q is not equivalent to the real 
world probability measure P. 



In a complete market, when assuming that Ft is the numeraire portfoho, 
one has for the corresponding Radon- Nikodym derivative A^ at time t the 
expression 

A? = (£) ', (6.6) 

which forms under the above 3/2 volatihty model the inverse of a squared 
Bessel process of dimension 5 = — > 2 under P. It is well-known that for 
(5 = 4 this process is a local martingale but not a true {A^ P)-martingale, see 
Revuz and Yor (1999). More precisely, A*^ is in this case a non-negative 
strict {A^ P)-local martingale, and thus, a strict {A, P)-supermartingale. 
The above 3/2 volatility model is then an example of a viable financial 
market model, as discussed in Loewenstein and Willard (2000), where the 
traditional notion of no-arbitrage, see Delbaen and Schachermayer (1998), 
cannot be verified. Therefore, a more general pricing method than the 
classical risk neutral one is needed to price derivatives. 



7 Real World Pricing 

It is now the aim to price derivatives under the real world probability 
measure. Since an equivalent risk neutral probability measure does not 
exist under the above model, one can follow the ideas in Platen and Heath 
(2006) and use the numeraire portfolio as numeraire or benchmark. For 
this purpose, assume that the diversified index Ft represents the numeraire 
portfolio and the savings account is the only other traded security to sim- 
plify matters. The numeraire portfolio is defined as the portfolio that when 
used as benchmark makes all nonnegative benchmarked portfolios supper 
martingales. This means, the benchmarked savings account is in our setting 
the inverse of Ft, and one has by (6.1) 



4 = :^(''-"'""*-iT')*"^'"'" <'■" 

It is well-known that as long as the drift in (7.1) is not strictly positive -^ 
forms a super-martingale. Therefore, let us request that 

r]>h (7.2) 

and 

ai > 27. (7.3) 

We will see below that these parameter constraints are very realistic. 

In this context the following notion turns out to be crucial: 



Definition 7.1 A price process U = {Ut,t G [0,T]}, with E 



\Ut\ 
Ft 



< 00, 



is called fair if the corresponding benchmarked price process U = {Ut 

Ft- 



%-,t G [0,T]} forms an {A, P) -martingale, that is, 



Ut = E 



Ut\A 



T\-^t 



(7.4) 



for alio <t<T < T. 



As discussed in Platen and Heath (2006), the minimal supermartingale, 
which replicates a given benchmarked contingent claim, is the correspond- 
ing martingale. Since the minimal price is economically the reasonable 
price for a replicable claim, it is the fair price that should determine the 
value of a derivative if no other constraints exist. For a replicable contin- 
gent claim Hf, payable at time T G [0,T] with E{'-p^) < oo, this yields 
the real world pricing formula 



T 



Vt = FtE(^^\A?j (7.5) 

for alltG [0,T],Tg [0,T]. 

Now, one can discuss the link between real world pricing and classical risk 
neutral pricing. As shown in Platen and Heath (2006), and as already 
indicated in (6.6), in a complete market when the Radon-Nikodym deriva- 
tive process A*^ = {A^,t G [0,T]} for the putative risk-neutral probability 
measure Q is given by the ratio 

At -^U- jr, (7.6) 

then, one obtains from the real world pricing formula (7.4) the equivalent 
expression 

yt = E(^^Hf\A?j (7.7) 

for all t G [0,T],T G [0,T]. Recall that zero interest rates are assumed, for 
simplicity. 

In the described special situation of a complete market the Radon-Nikodym 
derivative process A'^ equals the normalized benchmarked savings account. 
If the savings account would be a fair price process, that is, a martingale, 
then the candidate Radon-Nikodym derivative process A'^ would be an 
[A, P)-martingale. This would guarantee that the risk-neutral probability 
measure Q exists. In this case, one would obtain from (7.6) by the Bayes 
rule the standard risk neutral pricing formula with Vt equal to E^{Hf\At). 
However, when looking at discounted S&P 500 total return index data, 
one observes that A^ exhibits in the long run much smaller values than at 
the beginning. This is an obvious reflection of the existence of the equity 



premium. It suggests that it may not be realistic to model A*^ as an {A, P)- 
martingale when pricing derivatives over long periods of time. 

Note that in the above 3/2 volatility model with Ft as discounted numeraire 
portfolio, the Radon-Nikodym derivative is, in general, not a martingale 
because it is only a supermartingale. Moreover, it follows from (7.6) by the 
supermartingale property of the Radon-Nikodym derivative that 

A? > ^t(A^) (7.8) 

for < t < T < T. For a nonnegative contingent claim Hf, when re- 
expressing (7.7) by using (7.8) one obtains the inequality 

E(§Hf\At) 

Vt < \\o X • (7-9) 

t 



The right hand side of the above inequality could be interpreted as the 
formal " risk neutral" price, which can be substantially greater than the 
real world price. The concept of real- world pricing generalizes classical risk 
neutral pricing. It does not impose the restrictive condition that A*^ has to 
form an [A., P)-martingale. As a consequence, real world pricing removes 
from the assumptions of Carr and Sun (2007) and Itkin and Carr (2009) 

2 

the necessity to require the condition q < \, which is imposed by these 
authors to prevent the risk neutral volatility from exploding. The gained 
freedom allows us to focus on the modeling of the real world dynamics of 
the 3/2 volatility model. 



8 Analytical Formulas of Variance Swaps 

The evaluation of the real world price K;(t, Ft) of the variance swap for a 
discounted diversified equity index Ft at time t = is given by: 

K(0,Fo) = FoEl-"^ °'^ "'' 



= F,Le(^) - FoLK,E{^) . (8.1) 

Hence, this valuation of a variance swap can be reduced to the problem 
of calculating the expectation of the benchmarked underlying variance 

-2 

FqE{-^) and the fair zero coupon bond price Pt{0,Fo) = FqE{^). 

Since the contract is equal to zero at inception, the fair strike K^, or called 
the variance swap rate, is given by 

Ei^) 



For completeness, we shall present the transition density function 
Ps{(pt,Xt; ipf,Xf) of a time-transformed squared Bessel process 
X = {X^^, ift G [v^o, VfW of dimension 5 > 2, which refers to a move from 
Xt = X^^ at the transformed time ipt to the level Xf = X^^ at a later 
transformed time iff. From Revuz and Yor (1999), we have 



, , 1 lXt\'' I Xt+Xf \ ( ^Xt Xf 



2{Lpf-ipt) \xt J i 2(v9f-v9t)j Vv^f-V^t 

(8.3) 
for ift G [v?05 V'f ]) where I^ is the modified Bessel function of the first kind 



with index z/ = | — 1. Then for Lp G (0, oo) and 6 > 2 one can show that 



5 

the /3th moment 

OO for /3 < -f . 

Note that for /3 < — | the corresponding moment does not exist. The fair 
price of zero coupon bond Pt(^, Ft) at maturity T is 

Prit^Ft) = FtE{^\At) 



■F, 



T 



A;^Vexp{-6T}F^(r)''^V 
Ft „/ 1 



exp{-6T}AT^U^(T)'^* 

^* -(2(^(T)-^(t)))-iexpJ' "^^(*) 



exp{-6T}Ar^ ^"^ ^ "^ ^^^ ' \2{^{T) - ^{t)) 



7 



r(^) V7 ' 7 ' 2(v^(T) - v.(t)) 



^.5) 



where (fit) = f^(e — 1) and the function iFi(., ., .) is the confluent hyper- 
geometric function defined by 

Using (2.1), (3.1) and (3.6) the mean value of the benchmarked underlying 
variance can be rewritten as 

^'''-'\e(I11). (8.6) 



aT \ Y. 



Proposition 8.1 Let Y = {Yt : t > 0} satisfy the SDE 

dYt = (ai - bYt)dt + ^/2^tdW{t) (8.7) 

Let f3 = 1 + m - a + u/2, m = |(^ - 1), Yq = y > and u = 
-\/(a.i — 7)^ + 4/ii7. Then if m > a — | — 1, 



E 



(i 1 by — j^i^T f be 



-e 7(s''^-i) 



,6T 



^^!^^ ^^' r(/3) ^ . . . , ^ hy . 
72sinh2(f )y' r(l + uf^'^"' + ^' 7(e^^ - l)' 



Proof: Use Corollary 5.8 of Craddock and Lennox (2009) and observe that 
the fundamental solution in Corollary 5.8 reduces to the transition density 
of a square root process as ;Ui — )■ 0. Therefore, we have 

z-y{T,y,z)dz, (8.9) 

where 

h / r \ 

p{T,y,z) 



fz\^-< 2 



27sinh(6T/2) \y 

"~^ ~ tanh(6T/2) J 



f h f ^ , , 1/ + 2 \ 
X ^^P ( ^ I '^1-^ + (1/ - ^ ' — — ' 



''V7sinh(6T/2)/ 



(8.10) 



Then, by using formula 6.643.2 of Gradshteyn and Ryzhik (2000), if Re{m- 



q; + 1 + |) > 0, we have 



CO 



"1 1 



^ /-, , 1 ^\r f ^a/P 



V 27^ tanh(6T/2) 7 V7sinh(6T/2)/ 

r(/3) / b^ y' ( hy 
1 + z/ V27 sinh(6T/2) ; """^^ \2-i{e^T - V 

be'^ Y'""\, f by 



7(e''^-l); "-™-2'2^^(efeT_i)y' 

(8.111 



where Mfc^^(z) is the Whittaker functions of the first kind. 
Furthermore, from the fact that 

MkA^) = e"5^^5+M,Fi(i + ;, - A;, 1 + 2/i, z), (8.12) 

which is formula 13.1.32 of Abramowitz and Stegun (1972). Hence, along 
with some calculations, we have 

.g 7(e0-' -1) 



X 



2yym \^(e^'^-l)7^ 



U^sinh^(f)J f(rT^^^^(^'' + '^' 



7(e^^-l)'' 

(8.13) 



The result follows from the fact that 

"T ds 



E 



k 



Ys 



-^i?(e--^^""^17")u.=o. (8.14) 



D 



To further simplify the result of Proposition 8.1, the derivatives of the con- 
fluent hypergeometric function iFi(a2, &2 5 z) with respect to the parameters 
02 and 62 are summarized below, see Ancarani and Gasaneo (2008). De- 
noted G-i = G\a2M\z) = ^i^^^g^ and H^ = H^a^M.z) = ^iM|:^ 
the first derivatives of the confiuent hypergeometric function iFi{a2,b2,z) 
with respect to the parameters 02 and 62- 



G' 



_ Z^ V^ (Q2)mi(l)mi ^™"^ 

&2 ^ r62 + l)mi(2)™, mi! 

mi=0 / 1 \ y 1 

X 2F2{l,a2 + l+mi,2 + mi,b2 + l + mi;z) (8.15) 



Similarly, 



1 _ ^Za2_ sr^ (a2 + l)mi(&2)mi(l)mi ^""^ 
bl ^o(^2 + l)mi(&2 + l)mi(2)™imi! 



where (a2)„ ~ (^2+^) 



X 2i^2(l,a2 + l + mi,2 + mi,62 + l+"^i;-2) (8.16) 

t2+"-) 

r(a2) • 



The equation (8.8) can be rewritten in form of 
Jo Ys 



E 






bv 

1 , ,, ,-m+a-^u/2ru ,,/2-r(/3)r(l + Z>)^(l + Z>)/l4 



Mh2)-"'^''"'^\h 



20+lym ^V -' V V (r(l + Z>)) 



- -^h^(h2)-"'^''-'''Hhif'^ ,^^^^ (-G^ + H^)hi, (8.17) 

where <^ is a digamma function, h^ = ^a^inXi^) ' ^2 = (^i^h' ^3 = 

e"^^^^^'™^ /^4 = :^, /^5 = p^iFi(/3,l + z>,:^(^), ^ = f - 1, 
/3 = 1 + i> — a, G^ and H^ are given by (8.15) and (8.16) respectively with 
a2 = $,b2 = l + i' and z= ^^^^^ny- 



9 Hedge Ratio 

In this section, we shall discuss the hedging of variance swaps. Different 
methods on hedging variance swaps, have been proposed in the literature. 
These methods include the simple delta hedging, the delta-gamma hedging, 
hedging using option portfolios, hedging using a log contract and the vega 
hedging, etc. For a comprehensive overview of various hedging strategies, 
see Demeterfi et al. (1999), Howison et al. (2004). In Section 2, we 
show that the log contract does not work for the real world 3/2 volatility 
model. We shall derive the exact hedging formulas in this section. The 



delta hedging formula is derived in form of 



dE 



(T ds_ 

1 T-T 



dy 



y 



m+l 



{hif'^hM -m + 



yv 



by 



2/ii72sinh2(f) 7(e^^-l) 



by 



y™ 1 (1 + z/j 1 + z/ '-f[e°^ — 1) 

X (/i2)"™+""^/i4 (2 ln(2) + ln(/i2) - ln(/ii) - ^0) + ^(1 + z>) 

6^ 



+ ;^<"'>'''^'"{ - /^ 



T^sinh^lf 



-h4{h2)-'^+''-^/^h5 



2{h 



\-m+a-i>/2 ^ \P) i 



r 1 + Z> 



_62^_\mi 



L27(e^^ - 1)(1 + u) f^^ (2 + z>)^,(2)^, m,\ 



m,i=0 



6t/ 



X 2^^2(1,/? + 1 + "^1,2 + mi,2 + z> + mi, _^^^^ 



7(e^^ - 1) 



+ 



¥y 



-1)2 Z^ 



272(1 + v){e^^ - 1)2 f-^ (2 + z>)^,(2)„, mj 



mi=0 



X 2i^2(l,/3 + l + "^i,2 + mi,2 + z> + mi, 



61/ 



+ 



hy 



tti:. 



7(6"^ -1)' 
(/3)mi(l)mi (7(6"'^-!))'"' /3 + 1+mi 



27(l + z>)(e^^-l)^j2 + z>)„,(2)„, mj (2 + mi)(2 + z> + mi) 



X 2i^2(2,/3 + 2 + mi,3 + mi,3 + z> + mi, 



% 



7(6"^ - 1) 



hl5 



E 



^y wi 



+ l)mi(l + '>)mi)(l)mi (7(6"^-!)) 



(1 + ^)^(6^^-1)7^0 (2 + ^)„..(2 + z>)„.,(2)„, mi! 



X 2i^2(l,/3 + l + mi,2 + mi,2 + z> + mi, 
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E 



7(e^^-l)' 
+ 1)„,(1 + z>)™,(l)™, ^i(^(?^)™^"' 



(1 + ^)^(e^^ - 1)^7^ ^0 (2 + ^)-i (2 + ^)n., (2)™, 
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7(6"^ - 1) 



b(3y 
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mi! 



bg \mi 



(/3 + l)mi(l + «>)mi(l)mi (7(e6T'_l)) 
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(9.1) 



where z>= — — l,/3 = l + z> — a and 



/i 



7 



21n(2)/i4(/i2)"'"+"~"/'/i5 + (/i2)~™+"~'/'ln(/i2)/i4/i5 






In addition, the Vega of the variance swap can be calculated by 



dE 



where 



rT ds 
Jo Y7 



dE 



f 

JO 



T ds 



da 



is given by (9.1). 



dE 



(-v/27)y5. 



Jo 



ya 

1 rji 



dy 



(9.2) 



10 Numerical Results 



In this section, an example of variance swap rates is shown under the 
benchmark approach and the risk neutral approach. Recall the values of a 
parameter set in Section 4, we have b = rj, ai = b and 7 = 0.0137. 

Table 1 displays the prices of the variance swap rate of the benchmark 
approach and the risk-neutral approach for various maturities respectively 
and y = 0.3. 



Table 1: Variance swap rate under the benchmark approach 
and corresponding rates under the risk-neutral approach 



maturities 


rates under the benchmark approach 


rates under the risk-neutral approach 


0.25 


0.137767 


0.480415 


0.5 


0.18197 


0.434351 


1 


0.257293 


0.368659 


1.5 


0.31596 


0.323736 


2 


0.360989 


0.290785 



11 Conclusions 



This paper studies a 3/2 volatility model. It is shown that the risk neutral 
methodology can not handle this kind of 3/2 volatility model. However, 
the benchmark approach works well for this kind of model. Furthermore, 
the analytical formulas for pricing and hedging variance swaps are derived. 
In addition, the calibration and numerical example are demonstrated. 
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